Astronomy & Astrophysics manuscript no. 7326 


© ESO 2008 


February 5, 2008 





Letter to the Editor 



Accretion rate and the occurrence of multi-peaked X-ray bursts 



Anna L. Watts and Immanuel Maurer 



Max Planck Institut fur Astrophysik, Karl-Schwarzschild-Str. 1, 85741 Garching, Germany 



o 
o 



oo 



> 

m 
O 

o 

Oh 
i 

o 



ABSTRACT 

Most Type I X-ray bursts from accreting neutron stars have a lightcurve with a single peak, but there is a rare population of faint 
bursts that are double or even triple peaked. Suggested mechanisms include polar ignition with equatorial stalling, or multi-step 
energy release; the latter being caused by hydrodynamic instabilities or waiting points in the nuclear reaction sequence. We present 
an analysis of the accretion rate dependence of the multi-peak bursts, and discuss the consequences for the various models. The 
observations pose particular challenges for the polar ignition mechanism given current models of ignition latitude dependence. 
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1. Introduction 

Neutron stars in Low Mass X-ray binaries (LMXBs) accrete 
hydrogen (H) or helium (He) rich material from their compan- 
ion stars via Roche lobe overflow. This accumulating fuel burns 
and under certain circumstances the burning process becomes 
unstable, resulting in a thermonuclear flash known as a Type 
I X-ra y burst (for recent rev iews see IStrohmaver & Bildsten 
(2006); iGallowav et all d2007l) ). The nature of the burning de- 
pends primarily on accretion rate (|Fuiimoto. Hanawa & M ivaii. 
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Fushiki & Lambl [19871; iBildstenl 119981; iNaravan & Hevll 
Cooper & Naravaru [2006). At the lowest accretion rates, 
bursts are triggered by unstable H burning. Above ~ 1% of the 
Eddington accretion rate (M Edd ~ 1 X 10 _8 M o yr _1 ), H burn- 
ing stabilises. A shell of He builds up and this triggers pure He 
bursts. Above ~ 10% MEdd, not all of the H can burn before 
the He triggers, leading to various classes of mixed H/He bursts. 
Eventually, above ~ 30% MEdd, He burning also stabilises and 
no more bursts are seen. Observations support the existence of 
this upper limit dCornelisse et all 120031; IGallowav etaTll2007l) . 

The vast majority of X-ray bursts have a single peak in 
the lightcurve. However, there are some with a multi-peaked 
structure. These fall into two classes. The first are Photospheric 
Radius Expansion (PRE) bursts, where the flux reaches the 
Eddington luminosity. As the envelope is lifted from the star the 
temperature drops and the energy of the emitted photons falls be- 
low the X-ray band, leading to an apparent drop in countrate. As 
the envelope contracts again th e temperature rise s, and there is a 
second peak in X-ray emission dPaczvnskil [l983). If one looks at 
the bolometric rather than the X-ray luminosity, the double-peak 
structure vanishes. 

In this Letter we focus on the second class of multi- 
peaked bursts. These bursts have lower peak countrates 
and do not involve PRE, so that the multi-peak structure 
is also visible in the bolometric luminosity. Double-peaked 
bursts have been reported for 4 U 1608-52 dPenninx et al.L 
1989) i _G^17j^(|Kuulkers^LaU 120021) and 4U 1709-267 
dJonker et al.L 120041) . The largest sample, however, comes from 
4U 1636-536. Observations with EXOSAT revealed 5 double- 
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Fig. 1. 2-25 keV X-ray lightcurves for the four non-PRE multi- 
peaked bursts in the RXTE archive for 4U 1636-536. Burst 
numbers (l)-(4) are burst number s 40, 48, 56 and 111 respec- 
tively in the burst catalogue of Galloway et al.l d2007l) . PCU 
refers to the RXTE Proportional Coun t er Unit s. The burst mod- 
el led in | Bhattacharvva & Strohmaverl d2006al) is burst (3); that 
in Bhattachary ya & Strohmaverl d2006bl) is burst (4). 
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peaked bursts and one triple-peaked burst dSztaino et all [l985; 
Ivan Paradijs et all fl986; Le win et allll987l) . Subsequent obser- 
vations with the Rossi X-ray Timing Explorer (RXTE) have 
added four more d ouble-peaked bursts (Figure Q] see also 
IGallowav et all d2007l) ). 

Since the di scovery of these bursts, a number of cau ses have 
been suggested. iFisker. Thielemann & Wiescherl (|2004|) give an 
excellent review of the suggested processes and the observa- 
tional results that have led to some being ruled out. The most 
likely mechanisms that remain viable involve a stepped release 
of thermonuclear energy. This coul d be caused by hydrod ynam- 
ical instabilitie s leading to mixing (Fuiimot o et al.Lll988l) or, as 
suggested by Fisker, Thielemann & Wiescher (2004), a waiting 
point that impedes the nuclear reaction sequence. 
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More recent work has suggested that the critical factor may 
be the latitude at which the burst ignites. As will be discussed in 
more detail below, burst ignition is most likely to take pla ce on 
the stellar equator. Bhattacha rvva & Strohmaverl (2006a b) con- 
jectured that rare high latitude ignition could lead to a double- 
peaked burst if the burning front slows and partially stalls as 
it approaches the equator. This leads to a dip in the lightcurve: 
the second peak forms as the burning front recovers and accel- 
erates again. The authors illustrated the concept by modelling 
two of the double-peaked bursts of 4U 1636-536. A particularly 
attractive feature of this model is that it attempts to explain si- 
multaneously both the spectral evolution and the detectability of 
millisecond oscillations during the burst rise. 

So what determines ignition latitude? The general assump- 
tion is that accreting material comes into corotation and hy- 
drostatic equilibrium on a timescale short compared to the 
recurrence time between bursts. The material must therefore 
spread to form an equipotential surface that on a rapidly ro- 
tating neutron star is oblate, since the rotation acts to reduce 
the effective gravity at the equator. This means that the local 
accre tion rate varies with latitude, b e ing hi gher at the equa- 
tor. ISpifkovskv. Levin & Ushomirskvl d2002l) argued that this 
favoured equatorial ignition since it is here that the column depth 
necessary to trigger a burst is reached most quickly. 

The i ssue of ignition latitude h as now been revisited in more 
detail bv lCooper & Naravanl (120071) (hereafter CN07). These au- 
thors start by assessing the implic ations of a loca l accretion rate 
above which burning stabilises (Bildsten, 1998). As global ac- 
cretion rate rises, there must come a point where the equator has 
stabilised but higher latitudes have not. In this case, ignition will 
occur off-equator. The situation is complicated by the existence 
of different regimes of unstable burning, but detailed modelling 
by CN07 confirms the general principle. Ignition occurs at the 
equator until a certain critical accretion rate (the precise value of 
which depends on the stellar parameters and the accreted com- 
position). Above this rate, however, the ignition point moves 
rapidly towards the pole and remains there until the global ac- 
cretion rate is so high that all bursting activity is suppressed. 

The accretion rate dependence of the multi-peak bursts may 
therefore enable us to distinguish between the various mecha- 
nisms. If the multi-peaked bursts are indeed caused by polar igni- 
tion, for example, the results of CN07 indicate that we should ex- 
pect to find them at higher global accretion rates than the single- 
peaked bursts. In Section|2]we present an analysis of the accre- 
tion rate dependence of the multi-peak bursts. The implications 
for the various models are discussed in Section[3] 

2. Observations and Data Analysis 

To test the accretion rate dependence of the multi-peak bursts 
seen with RXTE we ca n make use of the RXTE burst catalogue 
dGallowav et al.L 120071) . hereafter G07. This catalogue records, 
for over 1000 bursts, a complete set of burst parameters includ- 
ing two that can be used to assess accretion rate. The first is po- 
sition in the colour-colour diagram. Sources like 4U 1636-536, 
know n as atoll sources, trace out a pat h in the colour-colour dia- 
gram (Hasinge r & van der Mi[l989). One can assign a colour- 
colour coordinate S z that traces progress though the diagram and 
is thought to be proportional to global accretion rate (Figure 4 of 
G07fl A precise calibration of this relationship, however, is dif- 

1 For atoll sources the colour-colour coordinate is usually called S a 
rather than S z (which is used for the higher accretion rate Z sources). In 
this paper, however, we use 5, for consistency with G07. 
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Fig. 2. Dimensionless persistent flux (as a fraction of the 
Eddington rate) against colour-colour diagram coordinate S z for 
X-ray bursts from 4U 1636-536, data from the RXTE burst cat- 
alogue (G07). S z is assumed to be proportional to accretion rate. 
Single peaked bursts are shown as diamonds, and the four multi- 
peaked non-PRE bursts by squares. Bursts (l)-(3) from FigureQ] 
cluster together; Burst (4) sits alone. 



ficult. The second measure that is used to estimate accretion rate 
is the persistent luminosity. This is considered to be less reliable 
than S z , but its advantage is that it can be used to calibrate the 
accretion rate as a fraction of M^&d provided that the distance is 
known (from, for example, a PRE burst). For a full discussion of 
these measures we refer the reader to Section 2.1 of G07. 

We start by considering 4U 1636-536, the source with the 
largest sample of multi-peak bursts. The RXTE catalogue con- 
tains 123 bursts from this source. Using PRE bursts as standard 
candles, G07 estimate that bursts are seen for accretion rates in 
the range 2-13% Meh- Figure|2]shows the distribution of single- 
peaked and multi-peaked bursts in a plot of persistent flux (nor- 
malised by the Eddington flux) against S z . 

Figure [2] shows that the multi-peaked bursts do not appear 
at the highest inferred accretion rates. There are, even for the 
most extreme burst (burst number (4) in FigureQ]) single-peaked 
bursts at higher accretion rates. The EXOS AT data for 4U 1636- 
536 support these findings. In May 1984 lSztaino et alj d 19851) 
made a day-long observation during which the persistent flux de- 
clined by « 30 %. Eight bursts were detected, four of which were 
double-peak non PRE bursts. Although the bursts were seen 
when the flux was relatively high, two sin gle-peak bursts were 
also d etected at this time. Later analysis bv lPrins & van der Klisl 
(119971) showed that at this time the source was in the middle ba- 
nana state (2.0 < S z < 2.5), not the higher accretion rate upper 
ba nana state. The late r double an d triple-peak e d burs ts reported 
by Ivan Paradiis et alj dl986l) and iLewin et all d!987l) were also 
seen in a sequence of over 20 other bursts, at a time when the flux 
was rising. Although both were seen towards the higher end of 
the observed flux, there are other bursts in the sequence at higher 
fluxes that are single-peaked. Colour-colour diagram position is 
similar to the 1984 observations: by this me asure too there are 
single -peaked bursts at higher accretion rates dvan der Klis et all 
119901) . 

What about the other sources that have exhibited non-PRE 
double-peaked bursts? For 4U 1608-52 there is one multi-peak 
burst in the EXOSAT sample dPenninx et all Il989t) and one in 
the RXTE catalogue (Burst 24 from this source, G07 Figure 9). 
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Fig. 3. The distribution of multi-peaked bursts (shaded) with S z 
compared to other bursts for the three sources in the RXTE cat- 
alogue that have (a) non PRE multipeak bursts and (b) trace out 
a long enough path in the colour-colour diagram to calculate S 
These sources are 4U 1636-536, Aql X-l and 4U 1608-52. The 
lower panel shows the distribution compared to all other bursts, 
the upper panel the distribution compared to PRE bursts. 



iKuulkers etall (l2002h and iJonker et all (12004 discovered one 
multi-peak burst apiece during RXTE observations of GX 17+2 
and 4U 1709-267. Figure 9 of G07 also reveals previously unre- 
ported multi-peak non PRE bursts from Aql X-l (Bursts 7 and 
32) and EXO 0748-676 (Burst 82). Every single one of these 
sources has single-peaked bursts at higher accretion rates. 

Given this sample, what can we say about the accretion rate 
at which multi-peak bursts occur? Three of the sources in the 
RXTE catalogue with multi-peak bursts trace out a full enough 
path in the colour-colour diagram to allow computation of S z : 
4U 1636-536, 4U 1608-52 and Aql X-l (although for the latter 
S z can only be calculated for one of the two multi-peak bursts). 
Figure[3]compares the distribution of the multi-peak and single- 
peaked bursts with S z for these three sources. Five of the six 
bursts cluster around S z = 2.1, but the remaining burst (Burst 
4 from Figure [U occurs at S z * 2.5. The multi-peak non PRE 
bursts do seem to be associated with higher accretion rates: a 
two-sided Kolmogorov-Smirnov test gives only a 4% probabil- 
ity of them being drawn from the same underlying distribution as 
the set of all other bursts. The distribution of multi-peak bursts 
is however completely consistent with the distribution of PRE 
bursts, which are also associated with higher accretion rates. 
However, the sample of multi-peaked bursts is small, so we 
should treat these statistics with some caution. 

We can also consider the dependence on absolute accretion 
rate as inferred from the persistent flux for those sources with 
a reasonable estimate of distance. 4U 1636-536 exhibits bursts 
at accretion rates from 2-13% MecM- Three of the multipeak 
bursts (Bursts (l)-(3) in Figure [1} occur at accretion rates in the 
range 3.2 to 3.7% Mem, while the other (Burst (4) in Figure[T|i 
occurs at 8% MEdd- 4U 1608-52 has bursts for accretion rates in 
the range 0.5-10% MEdd, with the multipeak burst occurring at 
3% Mndd- Aql X-l shows bursts for accretion rates in the range 
0.1-1 1% MEdd- The two multipeak bursts are at 2.3% M^id and 
3.3% Mndd- EXO 0748-676 has only been seen to burst at much 
lower accretion rates, in the range 0.5 - 1.5% MEdd- The multi- 
peak burst occurs at 0.9% MEdd- For 4U 1709-267 no calibration 



of accretion rate is possible as there are no definite PRE bursts. 
The multipeak burst from GX 17+2 is particularly unusual since 
this source is a Z rather than an atoll source and accretes at a rate 
apparently in excess of MEdd- However, the morphology of this 
burst does differ from the other examples. 



3. Discussion 

The recent modelling by CN07 indicates that polar ignition 
should only take place above a certain accretion rate. If multi- 
peak bursts are triggered by polar i gnition, as suggested by 
Bhatt acharvva & Strohmaverl d2006alfa) . it follows that they 
should occur at higher accretion rates than the single-peaked 
bursts. This is clearly not the case. The bursts of 4U 1636- 
536 also pose one other problem for t he polar ignition model. 
Bhattacharvva & Strohmaver (2006a b) model bursts (3) and (4) 
in FigureQ] For burst (3) they infer polar ignition, while for burst 
(4) they infer intermediate latitude ignition. Within the model of 
CN07, this would require burst (4) to appear at lower accretion 
rates than burst (3). The data show the exact opposite: burst (4) 
has a higher apparent accretion rate than bursts (l)-(3). So where 
does this leave the polar ignition model? 



3.1. Are multiple peaks a unique signature of polar ignition? 

The polar ignition model could still be viable if the majority of 
polar ignition bursts are single-peaked. In this scenario all of 
the bursts of 4U 1636-536 for S z > 2 would be polar ignition 
bursts. If this were the case, there are several issues to be re- 
solved. Firstly, this would set the accretion rate for the transi- 
tion from equatorial to polar ignition at a few percent of M E dd 
(G07, and Figure |2). The modelling of CN07 suggests however 
that this limit should be higher: in an example with a neutron 
star spinning at 650 Hz, CN07 find a transition accretion rate of 
» 16% M Edd . The spin r ate of 4U 1636-536 is lower, at 580 Hz 
(IStrohmaver et all fl998). which would push the expected transi- 
tion rate up. This disparity would need to be resolved. However, 
as we have discussed above, precise calibration of S z and accre- 
tion rate is difficult, and there are other areas of burst physics 
where similar discrepancies arise. 

We must then consider the reason why a few, but not the vast 
majority, of the bursts are mul t i-peake d: the stalling mechanism. 
Bhatt acharvva & St rohmaver (2006a b), revisiting an idea orig- 
inally suggested by Szt aino et al.l (119851) . proposed that stalling 
is caused by the interaction o f the burning front with a spr eading 
layer of accreted material (Inogamov & S unvaevl 119991) . They 
suggest that stronger bursts may be able to disrupt this spread- 
ing layer and overcome the stalling mechanism, giving rise to 
single-peaked bursts. If this is the case, we would expect the flu- 
ence of the multi-peak bursts to be lower than that of other bursts 
in the same region of the colour-colour diagram. Figure|4]shows 
that although the multi-peak bursts do have low fluences, there 
are single-peaked bursts at similar accretion rates with fluences 
that are the same or even lower. For the polar ignition model to 
survive it would be necessary to explain these low fluence single- 
peaked bursts. One possibility might be that these lower fluence 
bursts are so weak that the burning front stalls completely at the 
equator and never reaches the second hemisphere at all. 
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Fig. 4. Burst fluence against accretion rate as measured by 
S z (data taken from G07). Diamonds: single peaked non-PRE 
bursts; Asterisks: PRE bursts; Squares: multi-peak non PRE 
bursts. Error bars are similar in size to the diamonds. 



3.3. Alternatives to the polar ignition model 

If the polar ignition model cannot be reconciled with the ob- 
servation s, an alternative is required. The nuclea r waiting point 
model of lFisker. Thielemann & Wiescherl ((2004) predicts multi- 
peak bursts at accretion rates of a few percent MEdd. in line 
with the observations. The existence of multiple waiting points 
in the rp process also offers a natural explanation for the triple 
peaked burst, something that is hard to explain in the po- 
lar ignition model. The alternative is t wo stage burning trig - 
gered by hydrodynamical in s tabilitie s dFuiimoto et al.L [1988). 
Bhattacha rvva & Strohmaverl (l2006al) criticised this model on 
the grounds that it was hard to see how an unburnt layer of fuel 
could be maintained for long enough for the mechan ism to work. 
However, a recent study by iGallowav et alj d2006l) of the PRE 
bursts of 4U 1636-536 suggested the presence of an additional 
source of stratification in this star: this may promote the survival 
of unburnt fuel layers. Theoretical effort is required, however, to 
generate models of spectral evolution and timing behaviour for 
these alternative mechanisms in order to subject them to rigorous 
test. 



3.2. Do other factors affect the latitude dependence of 
ignition? 

The other possibility is that the models of CN07 do not fully 
characterize the problem. As the authors themselves note, whilst 
their models do reproduce, for example, average burst recurrence 
times, they do not necessarily reproduce the scatter around this 
average. There are cl early subtleties waiting to be resolved (see 
also the discussion in lCooper & Naravanl ( 120 06)). CN07 assume, 
for example, that variations in the conditions on the stellar sur- 
face are set only by the effective gravitational potential. We can 
break this assumption down into several elements. Firstly, that 
accreted fuel spreads rapidly and is not impeded in any way. 
Secondly that a similar re-adjustment occurs after a burst if any 
unburnt fuel remains. Thirdly, that there is no additional asym- 
metry that would affect burst triggering. If any of these condi- 
tions were violated, they might permit polar and equatorial igni- 
tion to co-exist at the same accretion rate. 

The assumption of f uel spread, based on estimates of its vis- 
cosity, is hard to dispute (Bild stenLll998l) . Magnetic confinement 
could lead to asymmetries in fuel distribution if the field were 
strong enough. However, there is no evidence for persistent pul- 
sations in 4U 1636-536, suggesting that the magnetic field is in 
this case too weak to channel or confine fuel. What about ad- 
ditional sources of latitude-dependent asymmetry? Two effects 
that should perhaps be considered are (a) the effect on igni- 
tion conditions (not just stalling condit i ons) o f a spreading layer 
at the equator dlnogamov & Sunvaevi 1 1 9991) . and (b) latitude- 
dependent temperature asymmetries. The latter could arise if, 
for example, burning fronts launched from the equator repeat- 
edly s talled when approaching the pole (due t o reduced flame 
speed. ISpifkovskv. Levin & Ushomirsky (2002)). Reduced heat 
injection into the crust at high latitudes could lead to a temper- 
ature gradient that might affect subsequent ignition conditions. 
Latitude dependent heating due to either core fluid oscillations 
or slightly asymmetric accretion could also play a role. We note 
that both of these mechanisms have been discussed as a means 
of generating a quadrupole moment that might limit the spin of 
rapidly rotating neutron stars via gravitational wave emission 
(Ushomirsky^taU|2000j). 



4. Conclusions 

We have presented an analysis of the accretion rate depen- 
dence of multi-peaked non PRE bursts, focusing in particular 
on whether they occur at higher accretion rates than single- 
peaked bursts. The fact that they do not contradicts one of the 
predictions of the polar ignition model, given the recent study 
by CN07 that predicts a strong dependence of ignition latitude 
on accretion rate. While there are ways in which the polar ig- 
nition/stalling model could be reconciled with the data, alterna- 
tive models should also be pursued. Ultimately a combination 
of models may prove to be the answer, particularly given that 
the multi-peak bursts are not a homogeneous set. To see this one 
only has to consider the differences in morphology and apparent 
accretion rate for the four bursts in Figure Q] 
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